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Abstract 

Experimental  results  and  analysis  of  the  effects  of 
boundary  layer  tripping  on  dynamic  stall  of  a  tran¬ 
siently  pitching  airfoil  are  presented.  At  low  Mach 
numbers,  the  tripped  airfoil  exhibits  qualitative  sim¬ 
ilarity  with  the  behavior  of  the  untripped  airfoil. 
However,  the  local  supersonic  flow  at  Mach  numbers 
greater  than  0.3  is  significantly  modified  by  the  trip 
leading  to  vastly  different  shock/boundary  layer  in¬ 
teractions,  dynamic  stall  onset  and  vortex  formation 
angles.  The  formation  of  the  laminar  separation  bub¬ 
ble  is  found  to  have  a  favorable  influence  in  delaying 
dynamic  stall  on  the  untripped  airfoil  flow.  In  both 
Mach  number  regimes,  the  tripped  flow  actually  stalls 
at  slightly  lower  angles  of  attack.  Further  experimen¬ 
tation  with  three  trips  on  an  oscillating  airfoil  showed 
that  the  dynamic  stall  process  is  very  sensitive  to  the 
state  of  the  turbulence  in  the  boundary  layer.  This 
sensitivity  points  to  a  need  for  finer  turbulence  model¬ 
ing  techniques  for  use  in  dynamic  stall  computations. 

Nomenclature 

Cp  pressure  coefficient 

c  airfoil  chord 

Dg  grit  size 

K  constant  based  on  wind  tunnel 

turbulence  level 

M  free  stream  Mach  number 
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1.  Introduction 

Dynamic  stall  is  a  problem  of  great  interest  to 
aerodynamicists  since  it  offers  the  potential  for  en¬ 
hancing  the  maneuverability  and  agility  of  aircraft. 
This  potential  is  significantly  affected  by  the  onset 
of  compressibility  effects  in  flows  over  pitching  air¬ 
foils.  A  sizeable  data  base1,2  now  confirms  that  above 
a  free  stream  Mach  number  of  0.3,  compressibility 
dominates  the  flow.  Firstly,  compressibility  promotes 
premature  stall.  Additionally,  the  large  flow  acceler¬ 
ation  can  cause  locally  supersonic  flows  and  induce  a 
series  of  shocks2  3  int  roducing  other  flow  interactions. 
Available  laboratory  experiments  on  compressible  dy¬ 
namic  stall  have  been  conducted  at  low  Reynolds 
numbers  (300,000  to  1,000,000-  with  the  exception  of 
McCroskey  et  al4  and  Lorber  and  Carta5).  These  low 
Reynolds  number  experiments  have  revealed  that  dy¬ 
namic  stall  originates  rapidly  over  a  very  small  angle 
of  attack  range3,6  as  the  laminar  separation  bubble 
that  forms  over  the  airfoil  bursts.  Further,  most  of 
the  events  of  dynamic  stall  are  concentrated  near  the 
leading  edge  region  of  the  airfoil  executing  the  rapid 
pitch-up  motion.  The  formation  of  the  separation 
bubble  clearly  indicates  that  the  flow  was  laminar  ini¬ 
tially  and  then  transitioned  to  reattach  as  a  turbulent 
shear  layer.  Flow  reattachment  generally  occurs3  be¬ 
tween  x/c  =  0.05  -  0.08  for  the  angle  of  attack  range  (6 
-13  degrees)  of  interest  in  dynamic  stall.  The  exact  lo¬ 
cation  depends  upon  where  transition  occurs  and  how 
the  transition  length  is  affected  by  the  local  flow  con¬ 
ditions,  in  particular  by  the  adverse  pressure  gradient. 
Gostelow  et  al'  point  out  that  the  transition  length 
is  reduced  considerably  by  large  adverse  pressure  gra¬ 
dients.  Thus,  it  should  be  expected  that  the  transi- 
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tion  process  varies  as  the  airfoil  pilches  up.  Since  low 
Reynolds  number  dynamic  stall  is  intricately  coupled 
to  the  bursting  of  the  separation  bubble,  the  process 
is  also  influenced  by  the  transition  process.  To  some 
extent,  even  high  Reynolds  number  studies  could  be 
affected  by  this  problem  since  the  boundary  layer  is 
laminar  near  the  stagnation  point.  Only  when  tran¬ 
sition  occurs  before  reaching  the  suction  peak  could 
one  expect  the  bubble  not  to  play  a  major  role. 

The  experimental  quantification  of  the  role  of 
transition  on  the  details  of  the  dynamic  stall  process 
is  a  daunting  task.  However,  it  is  also  an  important 
task  because  the  experimental  data  is  needed  for  es¬ 
timation  of  high  Reynolds  number  behavior  based  on 
low  Reynolds  number  experiments  as  well  as,  for  com¬ 
putational  code  validation  effort.  It  is  worthwhile  not¬ 
ing  here  that  owing  to  the  complex  nature  of  the  flow 
and  a  lack  of  the  understanding  of  the  physical  pro¬ 
cesses,  contemporary  computational  fluid  dynamicists 
either  ignore  transition  and  perform  a  fully  turbulent 
calculation8  or  arbitrarily  fix9  transition  at  a  certain 
x/c  location  and  invoke  a  turbulence  model  from  this 
location  in  the  computations.  Thus,  the  usefulness 
of  the  computational  results  becomes  severely  limited 
for  quantitative  comparisons.  Also,  extrapolation  of 
laboratory  results  to  prototype  or  flight  situations  is 
also  of  limited  value  until  this  issue  is  addressed.  In 
any  case,  validation  of  the  computed  data  needs  ex¬ 
perimental  results,  which  were  hitherto  not  available. 

The  present  paper  reports  some  experimental  re¬ 
sults  on  compressible  dynamic  stall  obtained  by  trip¬ 
ping  the  boundary  layer  with  a  three  dimensional 
roughness.  Most  of  the  results  will  address  the  behav¬ 
ior  of  dynam'c  stall  over  a  transiently  pitching  airfoil, 
which  were  obtained  with  a  roughness  height  of  about 
175pm.  Subsequently,  more  measurements  were  per¬ 
formed  to  investigate  the  flow  over  an  oscillating  air¬ 
foil  with  this  and  two  other  roughness  elements;  these 
will  be  compared  when  appropriate.  A  real  time  point 
diffraction  interferometry  (PDI)  technique  was  used 
to  obtain  the  global  as  well  as  surface  density  fields. 
The  images  were  processed  to  yield  pressure  distri¬ 
butions  to  compare  the  various  cases  quantitatively. 
The  results  clearly  demonstrate  the  need  for  proper 
tripping  of  the  boundary  layer  in  these  complicated 
flows. 


2.  Experimental  Facility  and  Techniques 


2.1.  Facility  Description 


optical  access  is  available  to  the  entire  flow  field.  Thus 
it  is  possible  to  study  the  flow  at  the  surface  near  the 
leading  edge,  where  the  dynamic  stall  vortex  forms, 
as  well  as  the  flow  field  away  from  the  airfoil. 

The  transient  pitching  motion  was  produced  by  a 
custom  designed  hydraulic  system.  Ref.  11  provides 
the  details  of  the  drive  and  these  are  briefly  given 
below. 


angle  of  attack,  a: 
pitch  rate,  d: 

maximum  acceleration  rate: 
change  in  a  during  acceleration: 
minimum  acceleration  time: 
free  stream  Mach  number: 
airfoil  chord: 

Reynolds  number: 


0-60° 

0-3600  °/sec 
600,000  0 f sec2 
<6°  of  pitch 
4  ms 
0. 1-0.5 
7.62cm 
2xl05  -  9x10s 


The  angle  of  attack  range  and  the  highest  pitch  rate 
of  3600  °/sec  on  the  7.62  cm  chord  airfoil  corresponds 
to  a  OO^/sec  pitch  rate  of  a  3m  chord  airplane  wing 
at  any  given  Mach  number;  thus,  the  rates  obtainable 
from  the  design  are  directly  applicable  to  flight  condi¬ 
tions.  The  system  controls  are  such  that  no  effects  of 
the  system  transients  on  flow  separation  are  notice¬ 
able.  The  system  uses  both  the  airfoil  position  and 
velocity  information  in  its  feed  back  loops  to  properly 
execute  any  pre-programmed  maneuver.  The  change 
in  angle  of  attack  during  acceleration  and  the  accel¬ 
eration  time  itself  are  limited  to  less  than  6°  and  4 
ms,  respectively. 

A  digital  optical  encoder  provided  the  instanta¬ 
neous  airfoil  position.  It  was  recorded  from  the  digi¬ 
tal  I/O  board  of  a  micro  VAX  II  workstation  and  the 
pitch  rate  was  timed  with  the  computer  internal  clock. 
The  pitch  rate  was  constant  to  within  1%  over  the 
angle  of  attack  range  of  interest  at  the  highest  pitch 
rates  since  it  was  specifically  tuned  for  such  rates.  At 
slower  rates,  a  mild  change  (less  than  2%)  was  noticed 
at  o  =  15°.  But,  since  this  angle  was  considerably 
higher  than  the  static  stall  or  the  dynamic  stall  onset 
angle,  it  is  believed  to  have  no  effects  on  the  results. 


2.2.  The  Trip 

Transition  was  fixed  by  applying  a  strip  of  ran¬ 
domly  placed  roughness  elements  along  the  upper  sur¬ 
face  of  the  airfoil,  near  the  leading  edge.  A  formula 
given  in  Ref.  13  was  used  to  estimate  the  minimum 
size  of  the  roughness  elements  required  to  trip  the 
boundary  layer.  The  formula  is 


The  results  to  be  presented  in  this  paper  are  part 
of  an  ongoing  dynamic  stall  research  program  in  the 
Navy-NASA  Joint  Institute  of  Aeronautics.  The  ex¬ 
periments  were  conducted  in  the  Compressible  Dy¬ 
namic  Stall  Facility  (CDSF)  in  the  Fluid  Mechanics 
Laboratory(FML)  of  NASA  Ames  Research  Center. 
The  facility  permits  study  of  dynamic  stall  over  a 
range  of  Mach  numbers,  using  non-intrusive  optical 
flow  diagnostic  techniques.  It  is  operated  as  a  part  of 
the  in-draft  tunnel  complex  at  the  FML  (for  details 
see  Carr  and  Chandrasekhara10).  In  the  CDSF,  the 
airfoil  is  supported  between  two  2.54cm  thick  opti¬ 
cal  quality  glass  windows  by  pins.  Since  the  pins  are 
smaller  than  the  local  airfoil  thickness,  unobstructed 


where  Dg  is  the  grit  size  in  inches,  R  is  the  Reynolds 
number  per  foot  of  length  for  the  free  stream  flow,  and 
K  is  a  constant  which  depends  on  tunnel  turbulence 
level  and  may  be  assumed  to  be  400.  This  formula 
indicated  a  grit  size  diameter  of  D,  =  56  -  89pm 
(0.0022  -  0.0035  in)  for  the  Mach  number  range  0.2  < 
M  <  0.3. 

Number  220  poIishCg  grit  with  a  size  range  of 
74  -  89  pm  (0.0029  -  0.U035  in)  was  selected  for  the 
construction  of  the  boundary-layer  trip.  The  grit  was 
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glued  to  the  airfoil  surface  with  a  water  soluble  adhe¬ 
sive  (Polaroid  print  coating  material).  The  region  of 
airfoil  surface  around  the  intended  trip  location  was 
masked  off  with  tape,  the  glue  was  brushed  on  to  the 
exposed  surface,  and  then  the  grit  was  applied  with  an 
air  brush  to  insure  an  even  distribution  of  the  rough¬ 
ness  elements.  The  trip  stretched  the  entire  span  of 
the  airfoil  and  extended  from  x/c  =  0.005  to  x/c  = 
0.03  along  the  upper  surface.  A  schematic  is  pre¬ 
sented  in  Fig.  1  indicating  the  location  and  thickness 
of  the  trip.  The  profile  was  obtained  from  a  digitized 
photograph  of  the  airfoil  under  no-flow  conditions  and 
indicates  that  the  maximum  thickness  of  the  trip  was 
175/im.  The  added  thickness  was  attributable  to  the 
adhesive  base. 

The  results  presented  in  this  paper  were  obtained 
for  a  transiently  pitching  airfoil  equipped  with  the 
trip  described  above.  Subsequent  experiments  were 
performed  on  the  flow  over  an  oscillating  airfoil  with 
this  and  two  other  trips.  The  oscillating  airfoil  data 
are  still  being  analyzed,  but  some  preliminary  results 
will  be  discussed  where  relevant.  A  spray-on  enamel 
lacquer  replaced  the  Polaroid  coating  as  the  bonding 
medium  for  the  grit.  The  lacquer  could  be  applied  in 
a  thiner  coat  and  had  a  longer  curing  time,  which  fa¬ 
cilitated  applying  the  grit  uniformly.  The  second  trip 
was  a  thin  version  of  the  first  trip  with  a  maximum 
thickness  of  ~  100/rm,  extending  from  x/c  =  0.005  to 
x/c  =  0.03.  The  third  trip  extended  from  x/c  =  0.05 
on  the  lower  surface  near  the  stagnation  point  and 
around  the  leading  edge  to  x/c  =  0.03  on  the  upper 
surface;  the  thickness  was  ~  130 pm. 

2.3.  Point  Diffraction  Interferometry  Tech¬ 
nique 

The  point  diffraction  interferometry  technique 
used  in  this  study  utilized  a  point  discontinuity  (in 
t  he  form  of  a  pin-hole)  located  at  the  image  of  a  point 
light  source  to  diffract  a  portion  of  the  incident  light 
into  a  spherical  reference  wave  front.  In  the  present 
application,  the  primary  optics  of  an  existing  schlieren 
system  were  used  (see  Ref.  12  for  details),  with  a 
pulsed  Nd:YACI  laser  replacing  the  conventional  spark 
as  the  light  source,  and  a  specially  created  point 
diffractor  replacing  the  usual  knife  edge.  The  laser 
light  was  expanded  through  a  microscope  objective  to 
fill  the  schlieren  mirror,  transmitted  through  the  test 
section,  and  refocused  by  another  schlieren  mirror. 
The  exposed  photographic  plate  used  to  create  the 
point-diffraction  spot  was  placed  at  the  focus  of  this 
second  mirror,  and  the  laser  was  pulsed  with  enough 
energy  to  burn  a  hole,  or  spot,  in  the  emulsion.  The 
spot  was  created  in  situ  by  passing  light  through  the 
test  section  at  under  no-flow  condition.  It  was  also 
precisely  tailored  to  the  application  under  investiga¬ 
tion,  automatically  correcting  for  nonuniformities  in 
the  light  source  or  optics.  The  tunnel  was  turned  on 
and  the  real-time  interference  fringes  were  recorded 
on  Polaroid  film(ASA  3000),  and  were  available  for 
immediate  viewing. 

The  analysis  of  the  interferograms  was  conducted 
using  a  package  developed  in-house  for  the  purpose. 
The  PDI  images  were  digitized  and  the  images  pro¬ 
cessed  on  an  IRIS  Work  Station  to  obtain  the  fringe 
interse-t  ons  with  the  airfoil  contour.  Using  isen- 


tropic  flow  relations,  the  fringe  numbers  and  hence, 
the  fluid  densities  were  converted  to  pressure  coeffi¬ 
cients.  This  assumption  was  made  even  through  the 
boundary  layer  and  the  vortex,  for  lack  of  better  flow 
relations.  However,  it  is  not  believed  to  introduce  sub¬ 
stantial  errors  in  the  pressure  field,  since  the  entropy 
change  is  generally  small,  until  deep  stall  occurs. 


2.4.  Experimental  Conditions 

The  experiments  were  conducted  on  a  7.62  cm 
chord,  NACA  0012  airfoil.  The  following  conditions 
were  chosen  for  the  tripped  airfoil  studies.  (The  cor¬ 
responding  data  for  the  untripped  airfoil  was  already 
available3). 


M 

a+ 

0 

0.02 

0.025 

0.03 

u.u3o 

0.04 

0.2 

X 

X 

X 

X 

0.3 

X 

X 

X 

X 

X 

X 

0.45 

X 

X 

X 

X 

The  chord  Reynolds  number  ranged  from  360,000 
at  M  =  0.2  to  810,000  at  M  =  0.45.  The  airfoil  was 
pitched  from  0  -  60°  about  the  point.  The  in¬ 
terferograms  were  obtained  at  the  desired  instanta¬ 
neous  angles  of  attack  during  separate  pitch-up  mo¬ 
tions  (one  picture/pitch-up)  for  0  -  20°  angle  of  attack 
by  strobing  the  laser  externally.  The  laser  pulse  was 
as  60  nanoseconds.  There  was  no  delay  between  the 
selected  angle  of  attack  and  the  angle  of  attack  at 
which  the  laser  pulse  was  actually  seen,  which  was 
detected  by  a  photodetector  that  latched  the  encoder 
display  when  the  laser  pulsed. 


2.5.  Experimental  Uncertainties 


The  following  are  the  estimated  uncertainties  in 
the  various  quantities; 


Mach  number: 
angle  of  attack: 
normalized  pitch  rate: 

Cp: 


jc, 

dir/c) 


±0.005 
0.1  degrees 

±0.5%  at  rates  >  2000°/.«fr 
±1%  at  rates  <  2000° /sec 
at  M  =  0.2,  ±  0.225 
at  M  =  0.3,  ±  0.075 
at  M  =  0.45,  ±  0.0375 
±10 


The  uncertainty  in  CF  is  estimated  to  be  1  fringe 
for  the  flow  in  general  with  about  3  fringes  possibly 
undetectable  for  the  peak  suction  pressure  coefficient. 


3.  Results  and  Discussion 


The  pressure  distributions  derived  from  an  analy¬ 
sis  of  the  interferogram  images  will  first  be  presented 
along  with  a  qualitative  discussion  of  representative 
images  used  as  flow  visualization.  A  discussion  of  the 
leading-edge  pressure  gradients  calculated  from  the 
pressure  distributions  will  follow.  Comparisons  will 
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he  made  between  results  obtained  on  an  untripped 
and  a  tripped  airfoil. 

3.1.  Interferogram  Images  and  Pressure  Dis¬ 
tributions 

Analysis  of  untripped  airfoil3  steady  -  and  un¬ 
steady  -  flow3  interferograms  of  the  airfoil  flow  field 
while  the  airfoil  rapidly  pitches  front  0  to  60  degrees 
angle  of  attack  has  shown  that  a  leading-edge  lami¬ 
nar  separation  bubble  forms  in  both  steady  and  un¬ 
steady  flows  for  all  Reynolds  numbers  of  the  experi¬ 
ment.  Furthermore,  for  the  low  Mach  numbers  stud¬ 
ied  in  the  untripped  case,  the  dynamic  stall  vortex 
originated  just  as  the  bubble  burst;  the  vortex  then 
grew  and  convected.  Tripping  the  boundary  layer 
ensured  that  the  flow  was  turbulent  from  the  origin 
of  the  trip,  which  means  that  a  laminar  separation 
bubble  could  not  form.  The  result  of  tripping  the 
boundary  layer  is  illustrated  in  Fig.  2  for  the  condi¬ 
tions  M  =  0.3,  o=ll°  and  a  nondimensional  pitch 
rate  of  a+  =  0.04.  Enlarged  interferogram  images  of 
the  leading-edge  region  for  the  untripped  and  tripped 
cases  are  presented  in  Figs.  2a  and  2b,  respectively. 
The  instantaneous  pressure  coefficient  distributions 
which  were  derived  from  these  images  are  given  in 
Fig.  2c.  For  both  cases  a  peak  suction  pressure  coef¬ 
ficient  of  Cp  =  —3.4  was  measured,  with  the  location 
of  the  suction  peak  in  the  tripped  case  shifted  slightly 
downstream  of  the  location  of  the  peak  in  the  un¬ 
tripped  case  (from  x/c  =  0.013  to  x/c  =  0.018).  The 
laminar  separation  bubble  appears  in  the  untripped 
airfoil  data  as  a  plateau  in  the  pressure  coefficient  plot 
which  extend;;  from  x/c  ss  0.025  to  x/c  ss  0.055. 

The  evolution  of  the  dynamic  stall  vortex  was 
greatly  affected  by  the  trip,  as  is  illustrated  in  Fig. 
3.  Again  for  the  conditions  M  =  0.3,  a+  =  0.04,  the 
earliest,  stage  of  dynamic  stall  vortex  development  be¬ 
gan  for  the  untripped  airfoil  as  the  airfoil  passed  in¬ 
stantaneously  through  the  angle  ct  =  15°,  as  can  be 
seen  in  the  interferogram  shown  in  Fig.  3a.  Here,  the 
early  stages  of  vortex  development  on  the  untripped 
airfoil  is  indicated  by  the  opening  of  the  separation 
bubble,  which  is  made  evident  by  the  change  in  cur¬ 
vature  of  the  fringes  at  the  downstream  end  of  the 
bubble  (compared  to  Fig.  2a,  where  the  fringes  of  the 
closed  bubble  meet  the  airfoil  normal  to  the  surface). 
In  contrast,  however,  as  the  tripped  airfoil  passed  this 
angle  (Fig.  3b)  the  vortex  had  already  grown  consid¬ 
erably  and  covered  the  first  25%  of  the  airfoil  chord. 
The  growth  of  the  vortex  can  also  be  seen  in  the  pres¬ 
sure  coefficient  distributions  shown  in  Fig.  3c.  As  the 
vortex  forms,  the  pressure  coefficient  distribution  be¬ 
gins  to  ‘lift,’  up  along  the  suction  surface,  and  as  the 
vortex  grows  in  extent  and  convects  downstream,  the 
increased  region  of  suction  along  the  upper  surface 
results  in  the  vortex  lift  of  dynamic  stall. 

The  development  of  the  pressure  coefficient  dis¬ 
tribution  over  the  angle  of  attack  range  of  a  = 
7°  to  a  =  17°  is  shown  in  Fig.  4  for  both  the  un¬ 
tripped  and  tripped  airfoils  for  the  conditions  M  = 
0.3  and  a+  =  0.03.  The  separation  bubble,  which 
began  to  develop  at  a  «  9°  on  the  untripped  air¬ 
foil  (Fig.  4a),  is  clearly  absent  on  the  tripped  air¬ 
foil  (Fig.  4b).  Further,  for  angles  of  13  degrees  and 


above  (Fig.  4c  and  4d)  the  untripped  airfoil  has  de¬ 
veloped  a  stronger  suction,  while  the  tripped  airfoil 
stalled  at  approximately  one  degree  lower  angle  of  at¬ 
tack  than  the  untripped  airfoil  (14°  vs  15°)  as  can  be 
seen  clearly  by  comparing  Figs.  4c  and  4d.  In  addi¬ 
tion,  the  location  of  the  suction  peak  moved  slightly 
downstream  for  the  tripped  airfoil  with  increasing  an¬ 
gle  of  attack,  unlike  that  for  the  untripped  case,  where 
it  moved  towards  the  leading  edge. 

At  higher  subsonic  free  stream  Mach  numbers 
the  leading-edge  flow  can  become  locally  supersonic. 
In  the  case  of  the  untripped  airfoil  it  was  shown2,3 
that  multiple  A— shocks  formed  near  the  leading  edge 
over  the  suction  surface  for  freestream  Mach  num¬ 
bers  above  M  ss  0.4.  In  contrast  however,  at  most 
only  two  shocks  were  observed  on  the  tripped  airfoil 
at  M  =  0.45,  indicating  that  forcing  transition  has 
dramatically  altered  the  boundary  layer/shock  inter¬ 
action  physics.  Interferograms  of  the  untripped  and 
tripped  flows  are  shown  in  Fig.  5a  and  5b,  respec¬ 
tively,  for  the  conditions  of  M  =  0.45,  a  =  1 1°  and  a 
nondimensional  pitch  rate  of  a+  =  0.025.  A  total  of 
five  shocks  formed  on  the  untripped  airfoil  for  these 
conditions  with  the  final  shock  occurring  at  nearly 
10%  of  the  chord  from  the  leading  edge.  Only  two 
shocks  occurred  in  the  tripped  flow,  with  the  final 
shock  located  near  the  downstream  end  of  the  trip  at 
about  the  3%  chord  point.  Fringe  counting  showed 
that  for  the  untripped  flow,  the  local  Mach  number 
was  supersonic  until  a  height  of  about  2%  chord  above 
the  airfoil  surface,  and  the  maximum  Mach  number 
was  about  1.25.  On  the  other  hand,  the  supersonic 
region  for  the  tripped  airfoil  extended  to  1%  chord 
above  the  airfoil  and  the  Mach  number  was  also  lower 
at  about  1.1.  In  both  cases,  the  supersonic  flow  tran¬ 
sitioned  to  the  outer  subsonic  flow  through  a  com¬ 
pression  wave,  through  which  the  interference  fringes 
were  found  to  bend.  The  Mach  number  where  this 
bending  effect  was  masked  was  about  0.85.  The  loca¬ 
tion  of  the  shocks  in  relation  to  the  trip  is  illustrated 
more  clearly  in  Fig.  5c,  a  composite  of  data  obtained 
from  digitized  images  of  the  trip  and  of  the  interfero¬ 
gram  shown  in  Fig.  5b.  The  bold  lines  in  this  figure 
represent  the  supersonic  flow  and  the  dashed  lines, 
the  remainder  of  the  compression  waves. 

Another  striking  feature  seen  in  both  the  inter¬ 
ferograms  of  Fig.  5  is  the  appearance  of  the  vertical 
fringes  following  the  last  shock.  In  Ref.  3,  this  was 
shown  to  be  associated  with  the  first  indication  of  the 
onset  of  dynamic  stall.  The  studies  reported  in  Ref.  3 
also  showed  that  even  after  the  dynamic  stall  vortex 
formed  (near  the  end  of  the  bubble),  the  leading  edge 
pressure  coefficient  remained  large,  until  the  influence 
of  the  vortex  was  felt  near  the  leading  edge.  After  this 
event,  the  pressure  coefficient  dropped  rapidly  as  the 
vortex  began  to  convect.  In  the  higher  Reynolds  num¬ 
ber  experiments  of  McAlister  et  a!14,  dynamic  stall 
onset  was  said  to  correlate  with  the  loss  of  the  suc¬ 
tion  peak.  However,  the  present  interferometry  stud¬ 
ies  seem  to  define  stall  onset  more  precisely.  Fig.  5b 
shows  the  presence  of  the  vertical  fringes  clearly  and 
a  slight  reversed  flow3  is  seen  for  0.03  <  x/c  <  0.06. 
Fringes  pointing  upstream  in  a  thin  layer  of  fluid  be¬ 
low  the  separated  shear  layer  are  evidence  of  the  lo¬ 
cally  reversed  flow.  This  leads  one  to  infer  that  that 
the  dynamic  stall  process  is  under  way  -  even  though 
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there  is  a  shock  still  present  above  the  shear  layer. 
On  the  other  hand,  although  the  vertical  fringes  are 
seen  for  i  lie  same  experimental  conditions  of  M  =  0.45 
and  iv  =  1  ln  for  the  untripped  flow,  dynamic  stall  has 
just,  been  init  iated,  since  the  region  of  reversed  flow  is 
still  not  easily  discernible.  Because  the  dynamic  stall 
vortex  originates  at  x/c  «  0.1,  a  definite  delay  ex¬ 
ists  between  when  the  vortex  forms  (either  because  of 
the  bursting  of  the  laminar  separation  bubble  or  due 
to  shock  induced  separation)  and  when  its  influence 
reaches  the  leading  edge.  During  this  time,  the  airfoil 
could  continue  to  develop  suction  as  the  angle  of  at¬ 
tack  is  increased.  This  difference  in  the  dynamic  stall 
onset  location  between  the  two  flows  explains  why  the 
tripped  airfoil  flow  stalls  at  a  slightly  lower  angle  of 
attack  than  when  it  is  not  tripped. 

3.2.  The  Peak  Suction  Pressure 

As  mentioned  earlier,  the  suction  pressure  reaches 
a  maximum  value  just  prior  to  dynamic  stall  onset; 
during  the  period  while  the  dynamic  stall  vortex  is 
growing,  the  suction  remains  high  and  then  drops 
sharply  once  the  vortex  begins  to  convect.  The  peak 
suction  pressure  coefficient  reached  is  shown  in  Fig. 
6  for  different  nondimensional  pitch  rates  and  Mach 
numbers.  Plotted  for  comparison  on  the  same  graph 
are  the  results  obtained  previously3  for  the  untripped 
airfoil,  ln  both  cases,  the  peak  suction  pressure  co¬ 
efficient  decreased  dramatically  with  increasing  Mach 
number,  due  to  the  increasing  compressibility  effects 
on  the  flow.  ThisefTect  is  seen  in  Fig.  6  for  the  tripped 
airfoil  as  a  generally  monotonic  increase  in  the  peak 
suction  pressure  coefficient  with  increasing  nondimen- 
•sional  pitch  rate.  For  the  untripped  airfoil,  however, 
tin'  peak  suction  pressure  coefficient  tended  to  be  rel¬ 
atively  independent  of  nondimensional  pitch  rate  for 
c*+  >  0.03( >  0.02  for  M  =  0.45).  This  independence 
of  the  peak  CP  on  nondimensional  pitch  rate  (for  high 
pitch  rates)  in  the  unt  ripped  flow  was  presumably  due 
to  the  presence  of  the  separation  bubble;  the  change 
in  the  shape  of  the  enveloping  streamline  due  to  the 
presence  of  the  bubble  altered  the  inviscid/viscous  in¬ 
teraction,  and  eventually  limited  the  growth  of  the 
suct  ion  peak  so  that  increasing  the  pitch  rate  further 
had  no  effect  on  the  peak  suction  level  reached. 

The  peak  Cp  tended  to  be  less  for  the  tripped  air¬ 
foil  than  for  the  untripped  airfoil  for  the  same  exper¬ 
imental  parameters  at  Mach  numbers  >  0.3.  At  the 
low  Macii  number  of  0.2  the  difference  is  represented 
by  a  single  fringe  in  the  interferograms;  therefore, 
within  the  accuracy  of  the  image  analysis  the  peak 
( j,  value  was  unaffected  by  the  trip  for  this  rase.  For 
the  higher  Mach  number  flows,  where  compressibility 
effects  dominated,  the  difference  was  significant;  five 
or  more  fringes  for  M  =  0.3  which  corresponded  to 
ACV,  >  0.81.  and  four  or  more  fringes  for  M  =  0.45, 
or  A  Cp  >  0.28.  The  differences  in  the  behavior  seen 
for  the  different  Mach  numbers  are  typical  of  this  flow 
since  compressibility  effects  set  in  at  M  =  0.32. 

3.3.  The  Adverse  Pressure  Gradient 

The  nondimensional  adverse  pressure  gradient 
following  the  leading-edge  suction  peak  was  deter¬ 


mined  from  the  pressure  coefficient  distributions  ob¬ 
tained  from  the  analysis  of  the  interferogram  images 
by  performing  a  linear  least  squares  fit  to  the  first 
three  or  four  data  points  following  the  highest  pres¬ 
sure  coefficient  value  measured.  As  an  example,  Figs. 
7a  and  7b  show  the  leading-edge  instantaneous  CV,  dis¬ 
tributions  on  the  suction  surface  of  the  untripped  and 
tripped  airfoils,  respectively,  for  the  angles  a  =  12°, 
13°,  and  13.5°  at  M  =  0.3  and  o+  =  0.03.  Indicated 
on  each  curve  is  the  straight  line  segment  used  to  es¬ 
timate  the  value  of  the  pressure  gradient  •  No 

fit  was  performed  for  the  a  =  13.5°  data  of  the  un¬ 
tripped  case  since  only  two  fringes  were  detected  in 
the  adverse  pressure  gradient  upstream  of  the  sepa¬ 
ration  bubble;  often  at  higher  angles  no  fringes  are 
detectable  in  the  suction  peak.  For  the  tripped  airfoil 
flow  (Fig.  7b),  the  dramatic  change  in  the  pressure 
distribution  and  the  gradient  which  occurred  between 
the  angles  a  =  13°  and  a  =  13.5°  was  a  result  of  the 
formation  and  convection  of  the  dynamic  stall  vortex. 

The  flattening  of  the  distributions  at  the  top  must 
be  commented  on  since,  in  general,  the  pressure  coeffi¬ 
cient  distribution  is  a  smooth  and  continuous  fun*,  it '>r. 
with  a  well  defined  peak.  Inherent  to  the  interferome¬ 
try  technique  are  in-plane  deflections  of  the  light  rays 
due  to  the  steep  density  gradients  found  at  high  an¬ 
gles  of  attack  (for  example  at  a  =  12°)  which  have 
been  shown  to  deflect  the  image  of  the  beam  into  the 
airfoil  surface15,  resulting  in  a  reduction  of  the  num¬ 
ber  of  fringes  visible.  These  ‘missing  fringes’  cause  the 
flat  distribution  and  result  in  an  under  estimation  of 
the  value  of  Cp  at  the  suction  peak.  It  is  also  possible 
that  the  trip  physically  obscured  some  fringes  very 
close  to  the  surface. 

The  leading-edge  adverse  pressure  gradients  for 
the  tripped  and  untripped  airfoils  are  compared  in 
Fig.  8a  for  the  nondimensional  pitch  rate  of  a+  = 
0.03,  and  M  =  0.2.  Also  shown  for  comparison  in 
Fig.  8b  are  the  peak  suction  pressure  coefficients 
for  each  angle  of  attack.  Until  the  separation  bubble 
formed  (at  o  =s  10°),  the  magnitude  of  the  CP  gra¬ 
dient  was  similar  for  both  the  tripped  and  untripped 
airfoils.  After  the  bubble  formed  the  value  of  the  gra¬ 
dient  rose  sharply  in  the  untripped  case.  For  both 
cases  the  gradient  attained  a  maximum  value  about 
two  degrees  (of  angle  of  attack)  prior  to  the  onset  of 
dynamic  stall  and  remained  relatively  constant  until 
the  onset  of  stall.  As  the  dynamic  stall  vortex  be¬ 
gan  to  form,  the  magnitude  of  the  gradient  dropped 
sharply.  This  sudden  drop  in  the  magnitude  of  the 
adverse  pressure  gradient  coincided  with  the  initial 
development  of  the  dynamic  stall  vortex  and  the  ap¬ 
pearance  of  the  vertical  fringes  in  the  interferogram 
images.  The  drop  occurred  at  a  =  14°  for  the  tripped 
airfoil  and  n  =  15°  for  the  untripped  airfoil.  Two  de¬ 
grees  later  the  leading-edge  suction  pressure  dropped 
in  magnitude  (Fig.  8b)  and  the  dynamic  stall  vortex 
began  to  convect  downstream. 

This  general  trend,  where  the  magnitude  of  the 
gradient  stops  growing  approximately  ‘"'o  degrees 
prior  to  the  onset  of  stall  and  then  drops  sharply  as 
the  vortex  begins  to  form,  appears  to  be  independent 
of  both  Mach  number  (Fig.  9)  and  nondimensional 
pitch  rate  (Fig.  10)  until  the  local  flow  becomes  su¬ 
personic.  Figure  9  shows  the  pressure  gradient  as  a 
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function  of  angle  of  attack  for  the  Mach  numbers  M  = 
0.2,  0.3,  and  0.45  and  the  nondimensional  pitch  rate 
of  0.03,  while  Fig.  10  shows  the  results  for  several 
pitch  rates  at  M  =  0.3.  It  is  seen  in  these  two  fig¬ 
ures  that  the  magnitude  of  the  adverse  pressure  gra¬ 
dient  from  which  the  flow  can  recover  increases  with 
increasing  unsteadiness  and  decreases  with  increas¬ 
ing  Mach  number.  This  is  consistent  with  the  pitch 
rate  and  Mach  number  effects  on  the  development  of 
compressible  dynamic  stall  which  have  been  reported 
elsewhere''3.  The  lack  of  data  for  the  untripped  airfoil 
at  higher  angles  of  attack  for  M  >  0.3  was  a  result  of 
the  PD1  image  distortions  (i.e.,  the  ‘missing  fringes’) 
discussed  earlier  in  this  section.  The  character  of  the 
pressure  gradient  at  M  =  0.45  is  of  particular  inter¬ 
est.  Fig.  9a  shows  the  untripped  airfoil  results,  where 
the  maximum  t  experienced  a  decrease  with  an¬ 
gle  of  attack,  at  most,  this  remains  constant 

for  the  tripped  airfoil  (Fig.  9b).  This  result  clearly 
demonstrates  the  strong  influence  of  compressibility 
on  the  flow  development,  reflecting  the  impact  of  the 
multiple  shocks  that  have  been  observed  at  these  con¬ 
ditions. 

For  all  the  parameters  investigated,  when  the  lo¬ 
cal  flow  remained  subsonic,  the  untripped  airfoil  was 
found  to  be  capable  of  withstanding  a  stronger  ad¬ 
verse  pressure  gradient  before  dynamic  stall  onset 
than  the  tripped  airfoil.  Likewise,  the  peak  suction 
pressure  at  stall  was  higher  on  the  untripped  airfoil 
than  on  the  tripped  airfoil  (refer  to  Fig.  6).  It  ap¬ 
pears  that  the  formal  ion  of  the  separation  bubble  had 
a  beneficial  effect  on  stall  for  the  conditions  of  the 
experiments  Evidence  for  this  conclusion  is  based 
largely  on  t!,<  data  shown  in  Fig.  8a,  which  clearly 
shows  a  similar  development  in  the  pressure  gradi¬ 
ent  for  both  tripped  and  untripped  flows  at  M  = 
0.2.  Moreover,  the  appearance  of  the  vertical  fringes 
(which  signify  the  onset  of  stall)  coincides  with  the 
drop  in  the  magnitude  of  the  pressure  gradient  where 
data  were  available,  and  occurs  two  degrees  before  the 
drop  in  peak  Cp  (which  occurs  as  the  vortex  convects 
downstream)  for  the  cases  studied.  This  would  indi¬ 
cate  that  the  mechanism  of  the  dynamic  stall  vortex 
development  was  unaffected  by  the  separation  bubble 
at  low  Mach  numbers.  The  picture  is  considerably 
different  at  higher  Mach  numbers  as  discussed  in  re¬ 
lation  to  Fig.  5. 

3.4.  Discussion  of  the  Experimental  Observa¬ 
tions 

The  results  presented  in  this  paper  briefly  de¬ 
scribe  the  first,  effort,  in  assessing  the  role  of  transi¬ 
tion  on  the  compressibility  effects  on  dynamic  stall  of 
pitching  airfoils.  It  was  originally  presumed  that  use 
of  a  transition  trip  would  render  the  local  boundary 
layer  turbulent  and  hence,  a  dramatically  different 
sequence  of  events  of  dynamic  stall  would  occur  and 
improved  stall  characteristics  for  the  airfoil  would  re¬ 
sult.  But,  contrary  to  the  expected  results,  a  detailed 
analysis  of  the  interferogram  images  revealed  that  the 
trip  for  which  the  data  are  presented  actually  caused 
dynamic  stall  at  lower  angles  of  attack  than  for  the 
untripped  airfoil. 


It  is  to  be  noted  that  as  shown  in  Fig  8a  at  M 
=  0.2  and  <v+  =  0.03,  the  pressure  gradient  when  the 
untripped  flow  first  stalls  and  forms  a  bubble  is  about 
100,  whereas  the  pressure  gradient  for  the  tripped  flow 
at  first  stall  (which  is  dynamic  stall)  is  about  200. 
These  levels  confirm  what  is  expected  of  the  ability  of 
the  tripped  flow  to  withstand  a  larger  adverse  pres¬ 
sure  gradient.  In  this  context,  the  value  of  350  for 
the  pressure  gradient  at  dynamic  stall  vortex  forma¬ 
tion  for  the  untripped  case  seems  apparently  incon¬ 
sistent.  However,  this  may  not  be  the  case,  since  the 
untripped  flow  originally  stalled  at  a  pressure  gradi¬ 
ent  of  100  forming  the  laminar  bubble  which  changed 
the  pressure  distribution  over  the  airfoil.  Dynamic 
stall  in  this  case  arose  out  of  the  ‘failure"  of  the  lami¬ 
nar  separation  to  reattach,  thus  bursting  the  bubble3. 
The  local  pressure  gradient  at  the  downstream  edge 
of  the  bubble  was  also  found  to  be  about  100.  Once 
the  stall  vortex  formed,  the  gradient  dropped  sharply. 
It  appears  that  after  the  bubble  forms,  the  pressure 
gradient  could  still  increase,  because  the  separated 
shear  layer  could  reattach  until  a  certain  higher  an¬ 
gle  of  attack.  The  pressure  gradient  at  stall  decreases 
with  increasing  Mach  number,  but  the  above  trend 
is  still  valid.  At  M  =  0.45,  the  vortex  forms  at  a 

very  low  value  of  j  =  25  in  the  tripped  flow. 
Data  could  not  be  processed  in  the  untripped  flow  be¬ 
cause  of  a  larger  region  of  the  supersonic  flow  and  the 
more  complicated  multiple  shocks/compression  waves 
interactions  affecting  the  process. 

Fig.  10a  and  10b  show  that  in  steady  flow,  the 
tripped  airfoil  stalls  at  a  lower  angle  of  attack  than 
the  untripped  airfoil.  (Similar  effects  have  been  ob¬ 
served  elsewhere,  cf:  the  results  for  the  NLR1  airfoil14 
where  the  tripped  airfoil  stall  is  1  degree  lower  than 
the  untripped  airfoil  stall  angle)  However,  even  for 
steady  flow,  the  values  of  the  adverse  pressure  gra¬ 
dient  when  the  laminar  separation  bubble  forms  and 
when  the  tripped  airfoil  stalls  are  very  close  (70-100). 
One  factor  that  could  affect  this  number  is  the  thick¬ 
ness  of  the  trip  used.  As  stated  earlier,  at  17.5/wj,  it 
was  considerably  more  than  the  estimated  boundary 
layer  thickness  of  about  60  /im.  The  trip  height  was 
based  on  recommendations  generally  given  for  a  flat 
plate  boundary  layer.  The  challenge  of  finding  the 
‘right  trip’  for  an  unsteady  flow  with  a  large  adverse 
pressure  gradient,  locally  supersonic  flow  and  in  which 
the  transition  point  moves  considerably,  is  formidable. 
Further,  standard  recommendations  on  trip  selection 
seem  to  be  based  on  the  fact  that  the  boundary  layer 
eventually  reaches  the  equilibrium  state.  When  lead¬ 
ing  edge  separation  occurs  as  in  the  present  dynamic 
stall  flow,  this  criterion  is  obviously  not  satisfied. 
Also,  the  stall  process  is  significantly  affected  by  the 
‘age’  of  the  turbulence  in  the  tripped  boundary  layer. 
If  the  turbulence  is  ‘young’  as  in  the  present  experi¬ 
ments,  major  differences  could  be  found  for  each  trip 
that  is  used.  In  fact,  the  experiments  were  repeated 
for  an  oscillating  airfoil  with  three  different  trips.  Pre¬ 
liminary  evaluation  of  the  data  showed  trends  similar 
to  those  described  in  this  paper,  although  the  pitch 
angle  was  varied  in  a  sinusoidal  manner.  Further  anal¬ 
ysis  is  on-going.  A  typical  plot  of  Cp  vs.  a  for  these 
trips  is  shown  in  Fig.  11,  in  which  changes  in  Cp  near 
stall  onset  caused  by  the  various  trips  are  clear.  This 
also  brings  out  the  extreme  sensitivity  of  the  dynamic 
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stall  process  to  the  local  flow  details.  This  conclusion 
will  be  of  great  significance  in  computational  model¬ 
ing  of  the  flow,  especially  when  selecting  a  turbulence 
model.  Until  the  computational  fluid  dynamics  mod¬ 
els  include  these  complex  physical  effects,  it  appears 
to  be  unreasonable  to  expect  quantitative  agreement 
between  experiments  and  computations. 

4.  Concluding  Remarks 

A  study  of  the  role  of  boundary  layer  tripping  on 
dynamic  stall  of  a  transiently  pitching  airfoil  has  been 
conducted.  The  following  conclusions  are  drawn  from 
the  study. 

1.  The  laminar  separation  bubble  present  in  the 
untripped  flow  was  found  to  have  a  beneficial  effect 
on  dynamic  stall  delay. 

2.  Dynamic  stall  onset  moves  closer  to  the  leading 
edge  in  the  presence  of  a  trip,  which  eventually  leads 
to  stall  at  lower  angles  of  attack  (by  about  1-1.5 
degrees)  than  observed  in  untripped  airfoil  dynamic 
stall. 

3.  The  leading-edge  adverse  pressure  gradient 
and  the  peak  suction  pressure  coefficient  were  lower 
in  value  on  the  tripped  airfoil. 

4.  The  behavior  of  the  flow  is  grossly  different 
under  compressibility  conditions  with  a  trip.  The 
shock/boundary  layer  interactions  are  modified  by 
the  trip,  as  also  is  the  leading  edge  pressure  gradi¬ 
ent. 

5.  The  sensitivity  of  the  flow  to  the  state  of  the 
boundary  layer  turbulence  points  to  a  need  for  highly 
refined  computational  flow  modeling. 
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Figure  5.  Comparison  of  Tripped  and  Untripped  Flows 
Over  the  First  15%  Chord,  M  =  0.45,  a+  —  0.025,  a  =  11°; 
(a)  Untripped  Flow,  (b)  Tripped  Flow,  (c)  Shock  Location 
Relative  to  the  Trip,  Corresponding  to  Fig.  (b). 


Nondimensional  pitch  rate,  a+ 

Figure  6.  Transiently  Pitching  Airfoil  Peak  Suction  Pres¬ 
sure  Coefficient  at  Stall  Vortex  Formation. 
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Figure  7.  Leading  Edge  Suction  Pressure  Distributions, 
M  =  0.3,  a+  =  0.03.  Slope  of  Solid  Lines  Indicates  Mag¬ 
nitude  of  Adverse  Pressure  Gradient;  (a)  Untripped,  (b) 
Tripped. 
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Figure  11.  Peak  Suction  Pressure  Coefficient  Development 
over  an  Oscillating  NACA  0012  Airfoil  for  Various  Bound¬ 
ary  Layer  Trips.  M  =  0.3,  k  =  a  =  10°  +  10°sm  2 xft. 
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Figure  10.  Effect  of  Pitch  Rate  on  the  Adverse  Pressure 
Gradient  Development,  M  =0.3;  (a)  Untripped  Flow,  (1>) 

Tripped  Flow. 


